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Abstract. It has been argued by Caron-Huot that infrared contributions to the jet quenching parameter in 
hot QCD, denoted by q, can be extracted from an analysis of a certain static-potential related observable 
within the dimensionally reduced effective field theory. Following this philosophy, the order of magnitude 
of a non-perturbative contribution to q from the colour-magnetic scale, g 2 T/iv y is estimated. The result is 
small; it is probably below the parametrically perturbative but in practice slowly convergent contributions 
from the colour-electric scale, whose all-orders resummation therefore remains an important challenge. 



1 Introduction 



One of the classic qualitative indications for the genera- 
tion of a thermal medium in heavy ion collision experi- 
ments is the fact that high-py jets get quenched [1,2]. In 
fact, not only do jets get quenched but they appear to 
do so more effectively than naive estimates suggest [3,4]. 
This motivates not only a complete leading-order weak- 
coupling analysis [5], but also developing methods to go 
beyond the leading order [6], taking into account large 
corrections from the infrared (IR) scales that are charac- 
teristic of thermal field theory [7,8]. (The weak-coupling 
regime has recently also been discussed in rcf. [9].) 

In a weakly coupled non-Abelian plasma, there are two 
momentum scales in the IR: the colour-magnetic scale, 
g 2 T/ir, and the colour-electric scale, gT (g 2 = 4ira s is 
the QCD gauge coupling). Wc refer to these through the 
gauge coupling, g E /ir, and the mass parameter, m E , of the 
dimensionally reduced effective theory [10, 11]. The contri- 
butions of the two scales have been disentangled for sev- 
eral thermodynamic observables, and a particularly rich 
source of information is the "spectrum" measured through 
screening masses [12]. For instance, the smallest screen- 
ing mass (at temperatures below ~ 10 GeV), which is 
parametrically of the form 2m E + 0(g E /ir), gets a numer- 
ically insignificant contribution from the colour-magnetic 
scale [13], whereas the Debye screening mass, which is 
parametrically of the form m E + 0(g E /n), is in practice 
all but dominated by the colour- magnetic scale [14]: 
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The purpose of the current study is to estimate whether 
one of these extreme scenarios could be relevant for q. 



2 Colour-electric contribution 

The rate of jet broadening (transverse momentum diffu- 
sion) is often parametrized by a phenomenological coeffi- 
cient, called the jet quenching parameter and denoted by 
q [15]. Although, strictly speaking, the definition of q is 
ambiguous even at leading order [16], it can be argued that 
the ambiguities can be hidden by considering a quantity 
called the transverse collision kernel, C(q ± ); then 
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It is now the choice of the upper bound q* which reflects 
the ambiguities. The form of C(q ± ) for q ± ~ ttT was de- 
termined in rcf. [5] and for q ± ~ gT in rcf. [6]; although 
the NLO contribution from q ± ~ gT to q is parametrically 
suppressed by 0{g), it is numerically large. 

More precisely, denoting by m E = gT^/N^/3 + N{/6 + 
0(g 3 T) the Debye mass parameter and by g 2 = g 2 T + 
0(g 4 T) the coupling constant of the dimensionally re- 
duced effective field theory, the next-to-leading order 
(NLO) expression for the IR part (q ± ~ m E ) of the colli- 
sion kernel can be written as [6] 
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Here C F = {N 2 - 1)/(2JV C ) is the Casimir of the fun- 
damental representation. Choosing m E <C q* <C 7rT, the 
integral in eq. (2) yields 

„ = glmlCp m q* gjm E C F N c f q* 
m : 2~ I 16m„ 
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The terms involving q* cancel against contributions from 
hard momenta, q ± > 7rT [5], and the q* -independent mid- 
dle term inside the curly brackets then represents the 
physical NLO contribution to q [6]. 



3 Colour-magnetic contribution 

For small momenta, q ± <C m E1 the NLO part of eq. (3) 
behaves as 



(5) 



We note from eqs. (3), (5) that, as also suggested by the 
operator definition [6], the small-<7 ± limit (q ± <C m E ) of 
C(<Z_i_) g° es over into minus the momentum-space static 
potential of three-dimensional pure Yang-Mills theory: 



C(q±) = -V(q ± ) , 



(6) 



where [17] 
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The Fourier transform to coordinate space, V(r) 



1)V(<Z ± ), yields 



V(r) = 
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where r* is a regularization-depcndent constant. However, 
at the next order the expansion breaks down: a direct 
momentum-space computation of V(q^) at 2-loop order 
produces in d = 3— 2e dimensions the (gauge- independent) 
expression [17] 
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A comparison with eq. (2) suggests the presence of a non- 
perturbative contribution to q at 0(g E ), both because the 
integral is logarithmically divergent at the lower end, and 
because the coefficient in eq. (9) is IR divergent. 



4 Transformation to configuration space 

In view the delicate nature of the IR contribution, it is 
useful to re-express the small-q^ part of eq. (2) in another 



way. The idea is to make use of eq. (6), in combination 
with a non-perturbative understanding of V(r), in order 
to get a handle on IR physics. 

To be concrete, let A (q±) be some cutoff function, 
where A is chosen to be formally in the range g E /^ <C 
A <C m E . Then we can rephrase the IR part of eq. (2), to 
be denoted by q A , as 

Q A = f^liC(q ± )e A (q ± ) = jd 2 vV 2 V{r)9 A {r), 

(10) 

where 8 A (r) = J (2^ e ~ tqj "' r ^a(i±) i s tne configuration 
space version of the cutoff function. For instance, if 
e A (q ± ) = e(A-q ± ), then 6 A (r) = AJ x {Ar)j2-Kr. This very 
choice is not particularly convenient, however, since the 
Bessel function J\ is oscillatory. We find it more transpar- 
ent to choose a Gaussian, 
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because both functions are elementary and positive. 

Now, because of rotational symmetry and the form of 
eq. (8), the behaviour of V 2 V(r) appearing in eq. (10) is 



V 2 V{r) = g 2 E C F d^(r) + 
= g 2 E C F S^(r) + 
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where we have adopted the notation F(r) = V'(r) and 
c(r) = r 3 V"(r)/2 from rcf. [18]. 

We subsequently need to subtract the known pertur- 
bative terms, eq. (8), from the IR contribution to q A , given 
that they were already included in eqs. (3), (4). Thereby 
we obtain an expression for the remaining IR contribution, 
let us call it 5q A : 
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where we rescaled the integration variable as f = rA and 
defined 



<P(r) 



F{r) | Mr) 7gjC F N c 



64tt 



(14) 



As eq. (13) shows, only the short-distance part of 4>(r), 
terms linear in r (modulo logarithms) , contributes to 5q A ; 
this corresponds to terms quadratic in r in V(r). 

The nature of the short-distance behaviour of V(r) can 
be discussed within a framework similar to the Operator 
Product Expansion [19,20]. The lowest-dimensional con- 
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densate in three-dimensional pure SU(3) evaluates to 
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where the errors come from Monte Carlo simulations 
(MC) [21] and a scheme conversion (NSPT) [22], and the 
numerical values apply for iV c = 3. Just cancelling the 
scale dependence from eq. (15) one might expect a short- 
distance behaviour of the form 



V(r) 
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However, according to ref. [23] the shape at g 2 r <C 1 is 
really more complicated, ~ g E r 2 ln[ln(l/g|r)]/ln(l/g , |r), 
where g 2 N c ln(l/g 2 r) represents the difference of octet 
and singlet potentials. Note that such a tail does not 
contribute in eq. (13) for A ^> g E /n. Unfortunately the 
analysis is not valid for the range ggr ~ 1 that is relevant 
for us here, so we resort to modelling in the following. 



5 Modelling of lattice data 

In fig. 1, numerical data for r 2 0(r) from ref. [18] is shown. 
The parameter ro is defined from 



rgF(r ) = 1.65 



(17) 



for any given lattice spacing [24]. In the continuum 
limit [18], 

glr « 2.2 , (18) 

and we have inserted this estimate in order to combine 
the last term of eq. (14) with lattice data. (Obviously, it 
would be nice to add data at shorter distances, but this 
task is non-trivial because strong cutoff effects appear and 
a careful extrapolation to the continuum limit is needed.) 

Now, motivated by the naive eq. (16) and the definition 
in eq. (14), we describe the data through the function 
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The data (at /3 = 20) are well modelled (x 2 /d.o.f. = 0.18) 
with a = 0.72(2), b = 0.55(1), c = 0.18(1). We remark that 
around r/ro « 2.0 this function is close to the asymptotic 
value given by the non-perturbative string tension minus 
the perturbative subtraction, {[0.553(1)] 2 — 7/16ir}g E r 2 ss 
0.81 [25]. In the following we assume that eq. (19) should 
reflect the magnitude of the true function for g 2 r ~ 1. 
Carrying out the integral in eq. (13), we get 
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Fig. 1. The function <fr defined by eq. (14), in units of the 
parameter ro, defined by eq. (17). The lattice data is from 
tables 5 and 7 of ref. [18]. 



If the ansatz of eq. (19) were correct, the coefficient a/rj], 
which comes together with a dependence on the cutoff A, 
should be an analytically computable function, cancelling 
against an NNLO contribution from the colour-electric 
scale. Within the model the coefficient b/r^ represents a 
"genuine" colour-magnetic contribution although, because 
of the logarithm, this notion is somewhat ambiguous. 



6 Phenomenological interpretation 

Omitting from eq. (20) terms vanishing for Atq 3> 1; re- 
placing the cutoff inside the logarithm with the scale ~ m E 
at which other physics sets in; inserting ro from eq. (18); 
and estimating m E /g 2 ~ 1 (cf. fig. 1 of ref. [13]); we get 
the order-of-magnitudc estimate 
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This can be compared with the middle term from the curly 
brackets in eq. (4), 
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The non-perturbative contribution, eq. (21), is clearly be- 
low the NLO perturbative contribution from the colour- 
electric scale, eq. (22). 

In conclusion, the contribution to q from the colour- 
magnetic scale, q ± ~ g 2 /7r, may well be of modest magni- 
tude. Thus the phenomenological motivation for its theo- 
retically consistent determination may be feeble; rather, 
it should probably be measured as a part of the total 
IR contribution from both the colour-electric and colour- 
magnetic scales, e.g. along the lines suggested in ref. [6]. 
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Note added 

Recently a paper appeared [26] in which, as acknowledged 
there, some basic ideas of the current study, communi- 
cated to one of the authors several years ago, were re- 
vealed. Unfortunately the practical implementation, citing 
e.g. a peculiar temperature dependence of Sq A , appears to 
suffer from misunderstandings. 
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